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Abstract 

The symmetry method is used to derive solutions of Einstein's equa- 
tions for fluid spheres using an isotropic metric and a velocity four vector 
that is non-comoving. Initially the Lie, classical approach is used to re- 
view and provide a connecting framework for many comoving and so shear 
free solutions. This provides the basis for the derivation of the classical 
point symmetries for the more general and mathematicaly less tractable 
description of Einstein's equations in the non-comoving frame. Although 
the range of symmetries is restrictive, existing and new symmetry solu- 
tions with non-zero shear are derived. The range is then extended using 
the non-classical direct symmetry approach of Clarkson and Kruskal and 
so additional new solutions with non-zero shear are also presented. The 
kinematics and pressure, energy density, mass function of these solutions 
are determined. 



PACs numbers: 0420, 0440, 0230, 0270 

1 Introduction 



Solutions of Einstein's field equations that describe the motion of spherically 
symmetric fluid spheres have been discussed by many authors for many reasons 
which include their wide ranging applications to cosmology and astrophysics. 
The most general solutions may be obtained using a metric in isotropic form 
for which the four velocity is non-comoving and which possess non-zero shear. 
However, such cases are in general far less mathematically tractable than their 
more restrictive comoving counterparts. One of the most important techniques 
available in the derivation of comoving, or shear free, solutions is the Lie, clas- 
sical point symmetry approach although the method has not been applied to 
non-comoving descriptions and non-classical or non-Lie approaches are seldom 
applied to the problem of solving Einstein's equations. 
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The Lie approach applied to comoving systems was originally described by 
Kustaanheimo and Qvist P and by Stephani 0, Stephani and Wolf [3] who 
present the full range of symmetries and show that these Lie solutions of Ein- 
stein's equations contain, two arbitrary functions of time. Of course there are 
many other known solutions involving only one arbitrary function of time per- 
haps most notably by McVittie @][5], but also by for example Nariai [S] and 
Chakravarty [7|- Whilst the mathematical links between the many solutions 
have been comprehensively discussed by for example McVittie 0, Srivastava 

Sussman and Stephani et al ^HI it is interesting to note that these have 
not been discussed within the overall framework of the symmetry method. Ul- 
timately, of course it is the physical properties of these solutions of considerable 
interest 

Solution's of Einstein's equations for fluid spheres described by a non- co- 
moving observer on the other hand are still relatively unknown inspite of recent 
interest by for example Ishak and Lake Ishak Senovila and Vera ^Sj 
and Davision [Tl|. Although such solutions often have kinematics rich in shear, 
acceleration and expansion the mathematical expression of the condition for 
isotropy of pressure is mathematically far less tractable than for a comoving 
system. Indeed the solutions of Narlikar & Moghe ^5]> McVittie and Wilt- 
shire mil discussed by Knutsen ^| ED E| were found with the aid of ad hoc 
mathematical assumptions. Clearly the symmetry approach provides a more 
structured environment in which to derive and discuss solutions of Einstein's 
equations with a non-comoving four velocity. 

It will be shown however that the complexity of Einstein's equations in the 
non-comoving frame reduces the range of Lie symmetries available compared 
with the comoving case. It follows that non-Lie approaches to symmetry such 
as the non-classical method of Bluman and Cole [201 and the direct method 
of Clarkson and Kruskal [22 should also be investigated to extend the range 
of similarity solutions.. It will be seen that the direct method based upon 
symmetries for comoving systems results in new solutions. 

Whist the focus here concerns the symmetry framework for solutions of 
Einstein's equations it is important to stress that the physical properties of 
all solutions must be established. For comoving cases these results are well 
known, see for example, McVittie [22] > Knutsen [221, Sussman |51 l24ll2"5] whilst 
Knutsen |17lll9j has examined some non-comoving solutions. A review of many 
physical properties has been discussed by Krasinski [2S] ■ 

This paper is organised in the following way. Einstein's equations, in terms of 
the isotropy condition for fluid spheres, are presented for the non-comoving and 
comoving four velocity in section 2. The kinematic parameters are evaluated in 
general terms. Section 3 describes the Lie symmetry approach, used in section 4 
to review symmetry solutions for comoving cases. The interrelationship between 
known solutions of Einstein's equations is also discussed. The Lie symmetry 
approach is extended to non-comoving cases in section 5 with a new solution 
presented. In section 6 non-Lie, non-comoving solutions with non-zero shear 
are derived using the direct approach of Clarkson and Kruskal and further new 
solutions, with their kinematic properties is given. 
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2 Fluid spheres, isotropy, kinematic parameters 

2.1 Preliminary equations 

It is the intention here to consider the isotropic coordinate system for which 

da 2 = e 2X dt 2 -e 2t - ' (dr 2 + r 2 dfl 2 ) dfl 2 = dO 2 + sin 2 (0) d<j, 2 (1) 

where A = X(r, f), /i = p,(r,t). In addition Einstein's field equations will be 
taken to be 

Gl = -8nTl Ti=(p + p) uW - Sip (2) 

where G\ is the Einstein tensor and TJ: is the energy momentum tensor of 
the fluid sphere. The pressure p, energy density p and mass function m,see for 
example, Misner and Sharp [52], Cahill and McVittie |2Hj may be calculated 
using: 

8np=G 2 2 %vp=G\-G\-G\ m = ^ {l + e 2 ^ _A VV? - (1 + rp r f] 

(3) 

where the sufhx r, t indicates a partial derivative. 

Furthermore u l are the components of the velocity four vector that satisfy 
u l Ui = 1 and where u 2 = = u 3 . These may be used to calculate the kinematic 
characteristics of the fluid JU]. In particular the expansion 9, acceleration v,i, 
shear tensor and shear invariant a are respectively: 

9 = u\ Ui = Ui-jU J 

Oik = U(i ; fe) - U(iUfc) - \h ik Q a = a lk a %k (4) 

where hik — gik — UiU k is the projection tensor and the parentheses indicate 
symmetrisation so that u^-k) = \ {u^k + Uk;i)- 

Using Walker (5^] the assumption of spherical symmetry for a perfect fluid 
leads to the isotropy condition: 

A - e 2X (Gjf + e 2 » (G 2 - G\) {G 2 - Gf) = (5) 

2.2 Non-comoving (non-zero shear) cases 

To determine the isotropy condition and kinematic parameters for noncomoving 
cases u 1 ^ it is convenient to write: 

G\-G\ = e-^Ki e 2X Gj = -e 2f "G\ = 2K 3 

G\-G\ = -2e- 2X K A - e- 2 ^K 2 (6) 

In this way the isotropy condition JSJ may be written as: 

A = e 2(f *- x) {AKl - 2K x K i ) - K Y K 2 = (7) 
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where if i, if 2, ^3 and if 4 are respectively: 

K\ = /x rr + X rr + X 2 — [i 2 — 2\ r fi r — — K3 = /i r t — \ r ^t 

r 

if 2 = /-J-rr — Kr — A r + /i r H K t = (J, tt — \ t fJ, t 

r 

From JHJ) and (0) and taking u 4 to be positive it may be shown that 

,1 = f = _^ll p 2X-2 t , 4 



(8) 



u = f = S e * " •« 1 7TT> (9) 



K 
2K 



(' 



2(A- M ) 



Thus after some calculation the expansion is 

9 = fX r + gXt + 2/ ^ rr + l) + / Mr + 3 5Mf + 9t + fr (10) 
r 

whilst the two non-zero components of acceleration are 

ui = -?y = -e 2 » (f\ r + / 2 /i. r + 2fg» t + f t g + ff r ) - \ r (11) 

The non-zero components of the shear tensor are: 

2g 2 e 2A+2 A1(/rAr + /rr _ /) 2fge^(f^ + f t ) 
au = 3~r 3 (12) 

Q-33 a 11 r 2 e- 2X 

a22 = ^Te = V - (13) 

O44 = (714 = 041 = (14) 

T 9 
and finally the shear invariant is 



(15) 



Equations (|12|) to (|15(l demonstrate that each component of the shear tensor 
is proportional to crnwhilst the shear invariant is proportional to a\ x . Hence it 
is necessary only to calculate a\\ to determine whether solutions of Einstein's 
equations have non-zero shear. 

2.3 Comoving (shear-free) cases 

In the case of a comoving observer u 1 = f = and u 4 = g = e~ A and so from 
HI 0|l to l|15|l the only non-zero kinematic parameters are: 

9 = 3/iie" A ui = -Xr (16) 
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and so the shear tensor is necessarily zero for metrics isotropic in form. For 
comoving systems the condition for isotropy of pressure is the particular solution 
of that gives rise to the two equations: 

Ai = K 3 = A 2 = K x = (17) 



3 Classical symmetry method 

Initially it is the intention to discuss solutions of these equations in the con- 
text of Lie's classical symmetry. In this method one-parameter infinitesimal 
point transformations, with group parameter e are applied to the dependent and 
independent variables (r, t, A,/i). In this case the transformations are 

r =r + ET] 1 (r,t, A,/x) + 0(e 2 ) t = t + erj 2 (r, t, A, fi) + O (e 2 ) (18) 
A = A + e^i (r, t, A, /z) + O (e 2 ) p, = fj, + e<j) 2 (r, t, X,fi)+0 (e 2 ) 

and the Lie method, see for example Stephani, 30 or Bluman and Kumei, |31| . 
requires form invariance of the solution set: 

£ = {A(r,t),M = /*(r,t), A = 0} (19) 

where A = {Ai, A 2 } is defined by l|17|) for a comoving observer and A = {A} 
defined by J7J in the non-comoving case. This results in a system of overde- 
termined, linear equations for the infinitesimals rji, rj 2 , </>i and cf>2 ■■ The corre- 
sponding Lie algebra of symmetries is the set of vector fields 

d d d d 

X =Vi K r , *> -V M) {r, t, A, jti) (x, t, A, /x) T^+'fe (r, t, A, ji) — (20) 

The condition for invariance of (|17f) and Q is the equation 

4 2) (A) | A=0 = (21) 

(2) 

where the second prolongation operator X E is written in the form 



4 2 )_A- + 4 t] — + n — ^, — ,,, — , — 

+ 4 ] ,f + 4^ + 4 t] 



5 . ,w 9 



- [**] 


d 






- 01 


dX tt + 


•01 




a 


+ 4 rt] 


a 













(22) 



where <^!*', ^-'',0- rr ',(/)- rt ' and 0-"' , i = 1, 2 are defined through the transforma- 
tions of the first and second partial derivatives of A and fi to the first order in 
s. Thus for example: 

Xf = X r + etffl (r, t, A, /i) jlf = Mr + e^jf' ( r i *> ^ A*) 

At = A t + e^* 1 (r, <, A, /x) fa = m + e(j}[ t] (r, t, A, /x) (23) 
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Once the infinitesimals are determined the symmetry variables may be found 
from condition for surface invariance applied to A = A (r, t), \i — /i (r, t): 

Oi=0i-7]iA r - r] 2 X t = ^2=02-?7lMr - V2Vt = (24) 

Both Maple and the now unsupported Macsyma software have been used to 
calculate the determining equations. 



4 Review of comoving (shear-free) solutions 
4.1 Symmetry reduction 

The determining equations resulting from the invariance condition l|21|) applied 
to (|17|) give the following solutions for the infinitessimals 771 ,772 , 0i and 02 in 
the generator 

Vi =Vi(r) = ^-+ c 2 r - -y- 0i = 0i (t) ( 25 ) 
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m = m (t) = i] (t) 02 = c + cir 

where Co, ci, C2 and C3 are constant symmetry parameters. Integrating the 
surface conditions gives 

A = +lnA(t) jti = $ (u;) + In (ft (r)) (26) 



-=/ r ln(ft(0)= / C °' C2+Cir2 3 rfr (27) 

J f + czr-2£ J V(t) J f +c 2 r-^f 

and A (t) is arbitrary. Using (|2(j|) , (|27|l and the functions z — z (uj), Z = Z (to) 
equation (|17|l is satisfied by the two ordinary differential equations 

z = e 9 = ^ = -^- (28) 

Zuu - zz^ - 2cqz w - z 3 - 2c z 2 + (2cic 3 + c\ - Cq) z = (29) 

and is the equation used by McVittie j^j- Thus his solutions are classical Lie 
solutions of the coupled system l|17|) . However i|29|) is also equivalent to: 

Z^ - 2c Z u + bZ 2 - (2cic 3 + c 2 - eg) Z = (30) 

where 6 is a constant. This equation may be generalised by referring to Kus- 
taanheimo and Qvist 1 who showed that solutions of Einstein's equations may 
be expressed in terms of L = L (x,t): 

L = e~^ x = r 2 e x = A(t)n t = -A(t)^- (31) 
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so that lfT7|l reduces to 

L xx = L 2 F(x) (32) 
Substituting the second of (|26[1 into l|32|l it follows that: 

F( X )=-b k -M - 2 =- h e~ C0l ^^ L (33 ) 

(cix 2 - 2c 2 x - 2c 3 ) ( Cl x 2 - 2c 2 x - 2c 3 ) 2 

and (|30|l is again recovered with Z defined by l|28|l . Equation (|33|l is a general 
form for F {x) and has been reported by Stephani and Wolf f3] in their symmetry 
analysis of i|32|) . Thus the Lie analysis of l|32|l admits solutions containing two 
arbitrary functions of time whilst solutions obtained by a symmetry analysis of 
the coupled system have only one arbitrary function, the similarity variable. The 
latter are the McVittie metrics which are otherwise identical. Hence Z = Z (uj) 
in (I30|l may be replaced by Z = Z (u>, t) which satisfies: 

Z^ - 2c Z u + bZ 2 - (2 Cl c 3 + c\ - of) Z = (34) 



so that 



e 



A _ AU\ Zt 



A(t)^- e» = L- 1 (x,t) = h(x)Z- 1 (u,t) (35) 



Z 



4.2 Relationships between u(x, t) and h (x) 

Equation ij2T|) gives rise to the particular cases that are presented in table 1 





h{x) 


LU (x, t) 


H{x) 


Para- 
meters Ci 


1 


e 2c 3 


x f dt 
2c 3 J TtfTj 




ci = 
c 2 = 


2 


Co - c 2 
H 2c 2 


r dt 
v(t) 


2 (c 2 x + c 3 ) 


ci = 
c 2 ^0 


3 


e c H 


tt r dt 


In ( c ^+^- C2 } 


ci ^0 
/?-c2 + 
2cic 3 > 


\Jc\x 2 — 2c 2 x — 2c 3 




4 


e c H 


J T](t) 


J 1— a I 

arctanl v-p J 


ci ^0 
£ = c|+ 
2cic 3 < 


V 1 c\x 2 — 2c 2 ir. — 2c 3 




5 


He Hc„ 




1 


ci ^0 
/3 = C 2 + 
2cic 3 = 


CiX - c 2 



Table 1: Relationship between h(x) and u> (x,t) with x = r 2 
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Note that when c\ = , cq = then from table 1 and l|33[) F (x) is constant 

giving the results of Wyman [35]. However when c\ = 0, cq = — 5c2/7 then 

_ 20 ...... 

F (x) = — 6 {— 2 [C2X + C3)} 7 giving the result of Wyman [231 • Also when 

co = 5c 2 /7, c 3 = 0, ci > then from ^ ( x ) = -bcj^x' 1 ^ ( Cl x - 2c 2 )~^ 
which gives the result of Srivastava [S] . 

4.3 Summary of solutions for Z(u,t) 

Consider now solutions of (|34[l required for the metric i|35[l . 

Firstly when 6 = then by (|33[1 F (x) =0 and i|32[) gives the simple form 
L (x, t) = A(t)x + B (t). Moreover (33 or J2HJ) with b = is wholly equivalent 
as has been discussed by McVittie [5] and Srivastava [Sj. 

In the case when 6^0 the solutions of ilML'l) are summarised in Table 2. This 
table expresses several solutions in terms of Weierstrass P-functions, P (y, t) = 
P (y, 92 (t) , 33 (*)) which satisfy 

p2 = 4F 3 - 52 F- 53 (36) 

where 32 (t), 33 (t) are the invariants of P(y,t). In the following C (t) is arbi- 
trary. 





Parameters 

Ci 


Z(w,t)= ff (a;)e-'' 


Weierstrass 
P-function 
properties 


1 


c = 


-§P( w ,t) + £ 


■92 = T2 

S3 = C(t) 


2 


c ^0 

^ 25 


256 9 ^ W ^ + 25b 


q [u)) = e 5 
52-0 

33 = C(t)^0 


3 


co^O 
= ^ 

25 


96c" (g+1) 

25b (g+2f 


g (w) = e 5 
52 = 0, g 3 = 


4 


CO 7^0 

,,2 

^ J 25 




g (o>) = e~ 
52 = 
ff 3 = C(t)^0 


5 


co^O 

„2 

5 = -2- 

l J 25 


24c* <f 
25b (g+lf 


g (w) = e s 
92 = 0,33 = 



Table 2: Solutions for Z (w, t) when b 7^ 
ruif/i /3 = 2cic 3 + c\ , x = r 2 



Using l|35|l table 2 (entry 1) gives rise to the metric 
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This is an explicit form of the solution also given by Srivastava |S] and presented 
implicitly by Stephani et. al. JUj .When C (t) is contant the solutions of table 
2 are those derived by McVittie 0. If g 3 = -/3 3 /216 then table 2 (entry 1) 
may be integrated to give the solutions of McVittie 0] . Solutions for the case 
93 = /3 3 /216 were discussed by Chakravarty et. al. [7]- 

The metric (|37|) may also be considered using the Jacobi elliptic function 
sn (u) = sn (u, k) with elliptic modulus k, < k 2 < 1. It may be shown that 

Z (to (x, t) , t) = ^ j- gn2 3 {u k) +l + k 2 ±Vk*-ki + 1 j (38) 

where h = h(t), the elliptic modulus k = k (t) and it is the independent variable: 

& ■ . g 44 (*• + !)(«"- 



"4VA: 4 -fc 2 + l 27 

(39) 

Solutions l|38|l with constant C (i) have been discussed by Sussman [?JJ. In 
particular when /3 > , fc = 1 then sn(it, 1) = tanh («) and solution of 0] is 
recovered whilst for k = then sn(u, 0) = sin (u) giving the solution [7]- 

Table 2 (entries 2, 4) apply for all values of cq ^ 0. Thus when C3 = then 
(3 = c\ and entry 2 gives co = ±5c2/7 discussed by Srivastava |Hj- When c\ = 
then again co = ±5c2/7 discussed by Wyman [33- table 2 (entries 2, 4) were 
first obtained by Wyman j^Sj directly from 1)32(1 and McVittie (HJ from equation 
l|29|) . Both applied techniques originally due to Ince Finally while entries 
2, 4 contain two arbitrary functions of time 77 (t) and C (t) table 2 (entries 3, 
5) contain only 77 (t) . Moreover the entries 2 to 5 are known only for particular 
values of j3 and so further investigation of the case Co ^ is necessary. 

5 Non-comoving (non-zero shear) solutions 

5.1 Lie symmetry reduction, general case 

Consider the symmetries of equation J7J using the condition for invariance (|21|) . 
In this case and without loss of generality the solutions of the determining 
equations may be written in the form: 

Vi = r m = V (t) (j>\ = c Q - T) t (j)2 = cq - 1 (40) 

The surface conditions (|24|) may now be integrated and it is found that 



A = * (lu) + cq j — - In 
J V 



rj (i = <f> (uj) + (cq — 1) In r 



w = lnr-|^ (41) 
Hence equation JJJ reduces to the ordinary differential equation 
e 2(*-* +C0W ) |( 2co ^ + K 5 f + K 6 (K 5 + 2c * w - 2c $ w )} - K 5 K 6 = (42) 
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where: 

K 5 = - V uu + $ 2 - * 2 + 2c $ w + Co - 1 

#6 - + - + - 2* w * u - 2c $ w - 2c tf w - c 2 + 1 = f^r 2 
rr)(t)K 3 = $ u y u -$ ww (43) 

Equation 142() was found orginally using ad hoc mathematical simplifications by 
McVittie and Wiltshire ^Hj- I n addition using J2J, the pressure, energy density 
and mass function satisfy 

8irp (r, t) = r- 2c ° {e~ 2 * [* UU) + * 2 + 

+e -2(*- CoW ) [ 2 $ U V U - 2$ ww - 3$ 2 - 2co$ w ] } (44) 

8tt P (r, t) = r- 2co je" 2 * [tf^ + * 2 - 2M U - 2c * w 

- 2$5 - 4c $„ - 2c 2 + 2] + 3 e -2(*-c^) $ 2 | (45) 

2m (r,t) = r c «e* | e 2 (*+ c o--*)$ 2 + 1 - ($ w + Co ) 2 | (46) 
Also with the additional substitution 

S (lu) = jf^rjy h (to) = v/4e 2 *+ 2c o"-5 2 (47) 

and from the non-zero components of the velocity four vector are 

U l = — ^-A- u 4 = ^(^ ( } 

Moreover the expansion is, u\ component of acceleration and an component of 
the shear tensor are respectively: 



e 



-<i> 



(H* u + 2H^ -H u + 2c Q H) S 



^ r co J{2 



r c 77 5 r c 5 

8 (2 if $0,-/^+ Co if) e -*+4*+4c ^ 2 $ w e-*+ 2 *+ 2c « w 



"' r//\S' 1 rS 

4e 2 *+ 2o ° w (if^ + if$ w -if w +coif) 

(50) 



rif 3 

32r co-2 e 5$+4 Co ^ (Hm^-H^) 8r c °- 2 e 3 ' I ' +2cotJ - $ w - c ) 



3H 4 S 3HS 
16r c «- 2 (if$ w - H w + c if) e -*+5<f+4c o, 

+ 3lf2 (51) 
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5.2 Some similarity solutions 



Solutions of Einstein's equation that exhibit a simple baratropic equation of 
state p = (7 — 1) p are of considerable interest and from (|44|l and (|45|l the 
resulting equation would be an ordinary differential equation in the variable ui. 
In principle one might attempt to solve this simultaneously with l|42[l . However 
it would not simultaneously be possible to fit the solution to empty spacetime, 
other than for dust. This follows because from l|44|) the condition j-'iol I36| for 
zero pressure at the boundary surface, r^ = r^ (f) implies that uj = uib where ujb 
is constant. However the boundary the mass function (|46|) should be constant . 
This is not possible except when cq = 0. Results corresponding to cq — and a 
particular solution of l|42|) for cq ^ are given in Table 3. 



Solutions of 
for $ and ^ 



Ancillary 
functions Q (uj) 
& conditions 



Constant 
parameters for 

E3> 



<f> = * - cqlu + In k 
InQH 



2a (k 2 - 1) 



Q = (A sinh y / au>+ 
BcoshVaw)e 2fc2acouJ 



c ^ 0, k ^ 1, a > 

20 = l+fc ^Lg-l ) 

Section \5.2.1\ 



$ = * - c q uj + In k 
lnQ(w) 



2a{k 2 - 1) 



Q = (A sin y/—au)+ 
B cos ^f^aui)e 2k2acouJ 



c ^ 0, k ^ 1, a < 

20 = l+fc '(2cg-l ) 

Section rntn 



$ = * - c w 
* = + Q (w) 



Q 



^uj + a a 



co ^ 0, k= 1 
Section UTFfl 



All $ and 4" 
^5-0 



T-solutions 



co-0 
ffecfion BH 



$ and \I> 

^5 _ 2(*-4>) _ 1 



R-solutions 

e 2(*-») _ 1 < 



c = 
Section [51HI 



Table 3: The functions <£> cmd ^ defining non-comoving, 
Lie solutions \4-l{ 

5.2.1 Kinematics for table 3 (entries 1,2,3) 

The non-zero components of the velocity four vector are: 

! e-^+^S 1 4 _ 2fc 77(f) e -*+ C0W 



where 



S = - 



2fe 2 e' 3/ (2c ^- 1) 
2 c fc 2 W w - 2 eg k 2 - 1 



= v^fc 2 e 2 * - S 2 



(52) 



(53) 



11 



In this way it may be shown that the kinematic parameters are: 

e t * tf -*(4fe 2 e 2 *S u -3f u S 3 +8)c 2 e"$ tf S) , 

e = n- — Tn L ( 54 ) 

£H A ui = 4fc 2 e 2 * (2* W S* 2 - 1SSJ) + 8fc 4 e 3 * (cqS-S u ) (55) 
+ 2fc 2 e* (*„- c ) S 3 -16fc 4 e 4 ** w 



16 fc 5 e°" 4 e 4 *" C0 w (5 + 2 e*) (S u -9 U S) 
<Jl ' 1 = 9^ (56) 



where: 



f wu = 2a (l-fc 2 )^ 2 +4coafe 2 ^ w - a^ 2 -^ (57) 

s 2fc 2 e* (4c 2 fc 2 * 3 -4c (c 2 fc 2 + l) fr 2 + (l-2c 2 fc 2 ) ^ + 2c 3 fc 2 +c ) 

(2c A: 2 ^-2c 2 A : 2 -l) 2 

(58) 

Note that solutions of table 3 (entries 1,3) have been presented by McVittie 
and Wiltshire with physical properties described by Knutsen ^| . The 
case of table 3 (entry 3) with A = gives rise to a stiff equation of state p = p 
and the metric may be transformed to a comoving system resulting in the static 
solution of Tolman |37|. Table 3 (entry 2) does not appear to have been been 
presented in the literature. 

5.2.2 Discussion of table 3 (entry 4), T-solutions, kinematics 
Consider H42fl with Co = with K$ = 0. Hence the metric is: 

da 2 = e™Tf 2 dt 2 - e 2 *r- 2 (dr 2 + rW) (59) 
From (gBJ) to (|4*5|l and applying K 5 = it follows that: 

* 2 (u 1 ) 2 ( U 4 ) 2 1 



2m(r,t) r 2 e -2$ 7y 2 e -2*-4* 



> (60) 



The expressions for l|44|l and l|45|l now depend only upon u>. On defining the 
timelike variable r and spacelike variable according to 

du) 

e 2(*~) _ l 



T = e * q = \nr+ I (61) 



then (|59|l may be transformed orthogonally to the form: 
da 2 = e 2V dr 2 - e 2W dq 2 - rW 



1 

$2 f e 2(«=») _ I) 



<' 2[ = 1 „., ... e 2W = e 2 * ( e 2 (*-*) - l) (62) 
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and so V and W depend on r alone. The isotropy condition G\ = G\ for the 
metric (|62|l : 

^ r+ ^- W -4 + ^^-^-=0 (63) 

The region of spaceime defined by l|62|) . (|63[) is the T-region, Novikov |38| . 
Some solutions of (|63|1 have been presented by McVittie and Wiltshire |39| . 

It may be shown using Q that the expansion, acceleration and non-zero 
components of the shear tensor for (|62[1 are given by: 

2 (tW t - 1) e 2W - y a 33 _ r (tW t - 1) e- v (64) 

CT 11 — n CT 22 — , 2 ~ — r 

3r sir 6< 3 

giving the shear invariant: 

2(rW T -l) 2 e^ 
o = ^ (65) 

Note that the vacuum solution of Einstein's equations corresponding to the 
metric (|59fl is given in terms of the Weierstrass P— function P (w) satisfying the 
equation P 2 = 4P 3 — + 36M 2 + -[g where 2M is the Schwarzchild radius. 
Using (|62ll the vacuum solution is thus: 

— = J-— = e 2VV = e- 2V = — - 1 (66) 

PI MM 2 (1-12P) 2 r W 

5.2.3 Discussion of table 3 (entry 5), R-solutions 

Now consider g2J> with c = 0. Hence K 6 - e 2 ^"*) (K 5 + K 6 ) = so that by 
O r/\i = -ZK^-qe 2 ^-^ . The metric is again (JSU but now (gBJ), O give: 



-2*^2 



_^ 2m (r, £) r 2 ?7 2 



> (67) 



On defining the spacelike variable, q and a timelike variable according to 

/duj 
! _ e 2(j>-*) ( 68 ) 

then the orthogonal metric is again (|59|) but now 

e 2V = ( 1 - e 2 (*-*)) e 2ly = — * TU (69) 

V / $ 2 (1 - e 2 (*^*)) 

and so V and depend on q alone. The isotropy condition G\ = G\ corre- 
sponding to the metric l|59l) is now the familiar 



1 (W q + V q ) 



e 



2W 



V qq + V 2 - V q W q - — - ^—2 y- + — = (70) 

q q z q g z 

The region of space-time defined by (52} ) <|70|) is the R- region, Novikov |38| . 
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6 Non-comoving, non-Lie cases with shear 



6.1 Introduction and symmetry reduction 

The non-classical method is a generalisation of the classical Lie group approach 
due to Bluman and Cole |2(J| that incorporates the surface invariant condition 
l|24|) into the condition of form invariance (|21[1 . Thus non-classical symmetries 
may be found by solving the non-linear set of determining equations 

4 2) (A) l A=0 ,o=o = (71) 

where $7 = {f2i,f22}-The non-classical method is one example of a more gen- 
eral conditional symmetry approach described by Ibragimov |4L)| in which the 
condition l|21|) for form invariance is supplemented by an additional condition. 

However, the complexity of this method and the condition of isotropy for 
an non-comoving observer Q) present prohibitive practical difficulties. Instead 
a restricted nonclassical approach, namely the direct method of Clarkson and 
Kruskal [2] will be employed. This lacks the generality of the non-classical 
method 01] but is easier to implement. The approach does not involve group 
theory but rather imposes a given similarity ansatz on the differential equations 
in question. In the context of this problem the ansatz i|26|) and 127|l for the 
comoving isotropy equations will be used. Thus: 

\ = -9{uj)+\nA(t) n = $ (to) + h (r) (72) 

where 

<*=[ * r-f-^r Hr)=f C °- C2 + Cir2 3 dr (73) 

will be applied directly to the general equation J2J. In this way the isotropy 
condition becomes: 

2K 7 tj 2 (t) r 2 (2$„*„ + 2co# w -2$ <J<J + irr-8c 1 C3) 

+ K 7 e- 2 *+ 2 *- 2 A+2 c ° / dt + 2 c « « @ u * u - - („ At + m ) 9 U ) 

+ 2 e- 2 *+ 2 *- 2 A+2 c " I TO dt+2 c ° " (*„ * w _ $_) 2 

+ 4c 3 K 7 r) 2 (* u -2$ w + 2cs-2cd) 

+ 2 Cl K 7 7i 2 r 4 (^-2$ a) -2c 2 -2co)=0 (74) 
where 

K 7 = 1^ + * 2 - 2 - 2 cqV u + (75) 

- $ 2 - 2 c $ w + 2 ci c 3 + c 2 . - c 2 Q 

Inspection of these equations reveals that r){t) ~ 1/t and rj(t) ^constant, give 
rise to solutions of Einstein's equations. These solutions are summarised in four 
cases presented in table 4. 
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6.2 Some particular solutions of (J74|) 

In table 4 c 3 = |, A = 1 are chosen without any loss in generality. 



2* 



,2<I> 



Ancillary 
functions 



Constant 
parameters 



kiQ 2 - c 1Q Q 



Q 



Hi -c 10 + c( H 



4c? 



= 



27?= J(c 5 w + c 6 ) 



-c 5 w 



- c 6 
1 

'2ci t 

cj = r 2 + c w t 2 



co,ci,c 2 = 0, 
c 3 = §, A= 1 
Section Unni 



2* 



2* 



»»(*) = 27 

= r 2 — t 2 



C0,Cl,C2 = 0, 
Section Wn\ 



e 2 "z 2 



z= (c 5 + c 6 e 2 ") 3 
e" = Vr 2 + c 3 e bt 



c ,ci = 
c 2 ,A=l, 
Section \6. 2.<H 



z 'i 



z -•• 



z = A + buj 

n{t) = -\ 



C0,Ci,C 2 = 

c 3 = \,A 
Section\fi 



= 0, 
= 1 

U 



TaWe ^: XTie functions $ and '5 defining non-comoving, 
non-Lie solutions of {7$ 

6.2.1 Kinematics for table 4 (entries 1 and 2) 

In the case of table 4 (entry 1) it may be shown that 

2c 5 (Hi -c 2 ) fli 
The pressure, density and mass function are found to be: 



87175 : 
87175 = 



Quj (5 /ci r 2 Q u - 


- 7c w ujQuj - 6ki Q 2 -\ 


- 6 cio Q) 




Q 3 (AiQ-cio) 




(5 fei r 2 Q Q u 


- 3 C10 W Qui + 6 fcl Q 2 


- 6cio Q) 




Q 3 (h Q - do) 




r 3 Qu, (ki r 2 QQ^ 


- cio u Quj + 2 fci Q 2 - 


2c 10 Q) 



2Qi (fci Q- cio) 
and from @ the non-zero components of the velocity four vector 

1 cio?7(i)u 4 4 r 

U = U = — ; = 

r yjQ (k x r 2 Q-c w uj) 



(76) 

(77) 

(78) 
(79) 

(80) 
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The expansion is 

2 c w t r (3 kl r 2 Q 2 - 2 ki do r 2 Q - Ak lCw uoQ + 3 c 2 w to) Q u 



e = 



3 1 

Qi {ki Q - cio) (fci r 2 Q - ci w) 2 



cio i (fci r 2 Q + 2 cio r 2 - 2 do J) 
ry/Q (ki r 2 Q - cio w) 5 
The component iii of acceleration is 

c w r (2£ 2 - w) (ki Q - cm) 

Mi = 2 

(fci r 2 Q - do w) 

r (2 fc 2 r 4 Q 2 - 3 fci ci w r 2 Q + 2 c 2 wr 2 - cf w 2 ) 

Q (fci r 2 Q - cio w) 2 

and the component an of the shear tensor is given through 



(81) 



(82) 



on = r (83) 

3 (fci r 2 Q - c w lo) 2 



where 



W = 2 fci cio r 4 Q w - 2 fci do w r 2 - 2 fc 2 r 2 Q 2 

+ 4 fci cio r 2 Q + fci cio w Q - 2 c 2 r 2 - c 2 10 oj (84) 

The case of table 4 (entry 2) with cio = — 1. is the conformally Minkowskian 
case where $ = ^ are arbitrary functions of to. However it has been shown see 
for example Infeld and Schild 02] an d Tauber 03| that the hyperbolic case 
fc = — 1 of the Robertson Walker metric may be transformed into table 4 (entry 
2) and so this case will not be pursued further. 

6.2.2 Kinematics for table 4 (entries 3 and 4) 

In the case of table 4 (entry 3) with C2 = 1 (without loss of generality) it is 
convenient to define the following additional functions. 

W = (6 r 2 7? + 3 c 3 z 3 - 4 c 6 e 2 w r 2 ) Y~ - 6 r 4 z 5 + 3 c 3 r 2 z 5 + 4 c 6 e 2 w r 4 z 2 
F = r 2 z 2 -6 2 e 2 "r 2 -c 3 6 2 e 2w (85) 

It may be shown that the pressure, density and mass function are 

4c 6 (9z 3 (r 2 + c 3 ) -5c 6 e 2 "r 2 ) 6 2 (27 z 6 + 36 c 6 e 2 w z 3 - 28 c| e 4 ") 



87TJ3 



9(r 2 +c 3 )z 8 9z 10 

(86) 
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b 2 (3z 3 + 2c 6 e 2 ") 2 4c 6 (9 r 2 z 3 + 9 c 3 z 3 + 5 c 6 e 2uJ r 2 ) 
87TP = 3^ 9(r 2 + c 3 ) z8 (8?) 



b 2 e 3uj r 3 (3z 3 + 2c 6 e 2w ) 2 2c 6 e 3uJ r 3 (3 r 2 z 3 + 3 c 3 z 3 + c 6 e 2uJ r 2 



18 (r 2 + c 3 ) 5 z 7 9 (r 2 +c 3 ) 5 z 5 

and from @ the non zero components of the velocity four vector are 



v 



4 



1 _b(r 2 + c 3 ) 



zW zy/Y 
The expansion is: 



(89) 



_b (l8r 2 z 3 Y + 9c 3 z 3 Y +I2c 6 e 2uj r 2 Y -W) 
3 r z 4 r 2 

while the non zero components of acceleration are 

(r 2 + c 3 ) o (-3 r 2 z 3 + 3c 3 z 3 + 10 c 6 e 2 " r 2 ) 
— — -Ui = 



Y 

(6r 4 z 4 -3c 3 r 2 z 4 ~4:C 6 e 2LU r 4 z) /2c 6 e 



Y 2 V 3z 

and the nonzero component an of the stress tensor is 



2u 

+ 1) (91) 



2be 2uJ rW , s 

o-n = (92) 

9(r 2 +c 3 )F5 

Finally the case of table 4 (entry 4) when C2= is a solution originally due 
to McVittie & Wiltshire ^B] with kinematics and other physical properties 
analysed by Bonnor and Knutsen |18j . 

7 Conclusion. 

It has been possible to show that many solutions of Einstein's equations for 
perfect fluid spheres are obtainable by the classical symmetry method. This is 
not only true for the well known spheres described by a comoving observer but 
it is also true for non-comoving cases. However these produce a more restricted 
range of symmetries and it proved necessary to employ non-Lie methods to ex- 
tend the range of similarity solutions. Application of the non-classical method in 
its full form proved impracticable for non-comoving cases. However a restricted 
implementation in the form of the direct method enabled certain new solutions 
of Einstein's equations to be found. The physical and kinematical properties 
of these solutions have been presented in most cases but a detailed analysis of 
their physical applicability will be considered elsewhere. 
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Whilst the classical symmetry method has been applied succesfully in the 
derivation of solutions for a comoving four velocity the analysis is far from 
complete as equation (|34|) . Abel's differential equation of the second kind re- 
quires further research. Moreover, non-classical techniques have not yet been 
applied to comoving systems and certainly more work is necessary into non-Lie 
approaches for non-comoving cases. In addition solutions (whether comoving 
or non-comoving) that admit a simple baratropic equation of state of the form 
p = (7 — 1) p are in general as elusive as ever. 

This paper has focussed upon symmetry reduction for line elements with 
spherically symmetry and certainly the methodology could be applied to find 
exact solutions for perfect fluids with for example, plane or hyperbolic symme- 
tries. These ideas will be discussed elsewhere. 
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